> 



X 



Radiative corrections to chiral separation effect in QED 

E. V. Gorbar/'^ V. A. Miransky,'^ I. A. Shovkovy,"* and Xinyang Wang^ 

^Department of Physics, Tarns Shevchenko National Kiev University, Kiev, 03022, Ukraine 
golyubov Institute for Theoretical Physics, Kiev, 03680, Ukraine 



'^Department of Applied Mathematics, Western University, London, Ontario N6A 5B7, Canada 

"^School of Letters and Sciences, Arizona State University, Mesa, Arizona 85212, USA 

''Department of Physics, Arizona State University, Tempe, Arizona 85287, USA 

(Dated: May 15, 2013) 

We calculate the leading radiative corrections to the axial current in the chiral separation effect 

in dense QED in a magnetic field. Contrary to the conventional wisdom, suggesting that the axial 

f-f\ , current should be exactly fixed by the chiral anomaly relation and is described by the topological 

contribution on the lowest Landau level in the free theory, we find in fact that the axial current 

f^ ' receives nontrivial radiative corrections. The direct calculations performed to the linear order in 

^v^ , the external magnetic field show that the nontrivial radiative corrections to the axial current are 

,_^^ ■ provided by the Fermi surface singularity in the fermion propagator at nonzero fermion density. 

^ ' PACS numbers: 12.39.Ki, 12.38.Mh, 21.65.Qr 

I. INTRODUCTION 

(— I ' Recently there was a significant interest in the dynamics of relativistic matter in a magnetic field. Assuming that 

Q_ii QCD topological fluctuations produce local V and CV-odd states \^ leading to a chiral chemical potential fj.5, it was 
' ' suggested that there exists a non-dissipative electric current j = e^B/i5/(27r^) in relativistic matter in a magnetic field 
br'l B [2431 ■ This phenomenon is known in the literature as the chiral magnetic effect (CME). (For a recent review see 
(— I . Ref. [a.) Moreover, the charge-dependent correlations and flow, observed in heavy-ions collisions at RHIC 0-01 and 
LHC [101 1 appear to be in a qualitative agreement with the predictions of the CME [ll|, [l2l ■ 



Unlike the chiral chemical potential, which is a rather exotic quantity and not so well defined theoretically, the 
chemical potential /z (associatcdj_for example, with conserved electric or baryon charges) is common in many physical 
\^ , systems. It was shown in Refs. J13l - [l5| that a non-dissipative axial current J5 = eB/i/(27r^) exists in the equilibrium 
f^ ' state of noninteracting masslcss fermion matter in a magnetic field. This effect is known as the chira l sep aration effect 
^^ [ (CSE) in the literature. (For a brief review, see Sec. 2 in Ref. 15|.) In fact, as suggested in Refs. [la, [131 j ^^^ CSE 
^T ■ may lead to a chiral charge separation (i.e., effectively inducing a nonzero chiral chemical potential /X5) and, thus, 
tH- ■ trigger the CME even in the absence of topological fluctuations in the initial state. 

f^ . The approach in Refs. [iJ: UM '^as based on the use of the operator form of the chiral anomaly relation [Tsj. It 

C^ '_ is well known that the corresponding relation, calculated at one-loop order, is exact and, as such, it cannot get any 
higher-order radiative corrections [T^. Therefore, it was argued in [ij, [l5| that like the chiral anomaly, the one- loop 
result for the axial current density, J5 = eB/i/(27r^), should be exact as well. 

Since the fermion propagator in a magnetic field depends nonlinearly on the magnetic field, the linear dependence 
of the axial current on B calls for a physical explanation. Using an expansion over the Landau levels, it was shown 
J^ I in Ref. [lj| that the axial current J5 = eB/i/(27r^) is topological in nature (see also Ref. [20| for a nice exposition 
and some details) and is defined by the fermion number density on the lowest Landau level (LLL). Moreover, it was 
shown |l4{ that a similar result holds even for massive fermions at finite temperature T, where the axial current equals 
js = eBni(TO, r)/(27r) and nL{m,T) is the effective one-dimensional (along the direction of magnetic field) fermion 
number density on the LLL. At zero temperature the axial current is given by J5 = eB-^///^ — m? /{2'k'^). Of course, 
in the chiral limit ttj — )■ this reduces to the same expression for the axial current as derived from the chiral anomaly. 
Note, however, that the connection between the induced axial current and the anomaly relation is not obvious beyond 
the chiral limit. 

The chiral anomaly is exact as an operator relation, but it contains the divergence of the axial current rather than 
the current itself. Consequently, to get the axial current from the chiral anomaly one should "integrate" the anomaly 
and calculate the ground state expectation value of the corresponding operator. Then, the question concerning an 
"integration constant" in the induced axial current and its dependence on interactions naturally arises. Until now, no 
conclusive answer to this question was given (e.g., see the discussion in Ref. ^). 

The first studies of the interactions effects were done in Refs. [la l2l| - |23| in the framework of the dense Nambu- 
Jona-Lasinio_(NJL) model in a magnetic field. Using the Schwinger-Dyson equation for the fermion propagator, it 
was found [la, [111, [221 that the four-fermion interactions generates a chiral shift parameter A. In the chiral limit, 
this parameter determines a relative shift of the momenta in the dispersion relations for opposite chirality fermions 



fc^ — )• fc'^ ± A, where the momentum k^ is directed along the magnetic field. The presence of the chiral shift parameter 
leads to an additional dynamical contribution in the axial current. Unlike the topological contribution in the axial 
current at the LLL, the dynamical one affects the fermions in all Landau levels, including those around the whole 
Fermi surface. Further, it was explicitly checked in Ref. [221 that although the axial current gets corrections due to 
the NJL interactions, the chiral anomaly does not. 

Since the NJL model is nonrenormalizable and the chiral anomaly is intimately connected with ultraviolet diver- 
gencies, in order to reach a solid conclusion about the presence or absence of higher-order radiative corrections to the 
axial current, one should consider them in a renormalizablc model. In the present paper, assuming that the magnetic 
field B is weak and using the expansion in powers of B up to linear order, the leading radiative corrections to the axial 
current in QED are calculated. We find that they do not vanish and attribute this result to the singularities in the 
fermion propagator at the Fermi surface. On the technical side, the iesign(fco) prescription in the fermion propagator, 
which is the only thing that distinguishes a chemical potential from the time component Aq of the photon field, plays 
a crucial role in deriving this result. 

This paper is organized as follows. In Sec. HI] we introduce the model and set up the notation. Also, we discuss 
some properties of the fermion propagator and the one-loop self-energy in the presence of an external magnetic field 
and a nonzero density. The calculation of the leading radiative corrections to the axial current is presented in Scc. lIIII 
We start from the formal definition of the current in terms of the fermion propagator, use its systematic expansion in 
powers of the magnetic field, and finally perform the explicit calculations. Our discussion of the results and conclusions 
are given in Sec. IIVI A new form of the Schwinger parametrization for the fermion propagator in the case of a nonzero 
magnetic field and a nonzero chemical potential, utilized in the main part of the paper, is presented in Appendix \X\ 
The details of the calculations of the radiative corrections to the axial current are given in Appendix [Bj 

II. FERMION SELF-ENERGY IN A MAGNETIC FIELD 

The Lagrangian density of QED in a magnetic field is given by 

C - -^F^-'F^, + 4, (i^V, + M7" - ™) '0 + '^2^(*7"5. + M7" + eA-;^'Y)4' ~ S^i^ij, (1) 

where /i is the fermion chemical potential, the last two terms are counterterms (we use the notation of Ref. [2J], but 
with the opposite sign of the electric charge, e — )- — e), and the covariant derivative is I?^ = d^ — ieA^ — ieA'r^*'. 
Without the loss of generality, we assume that the external magnetic field B points in the -I-X3 direction and is 
described by the vector potential in the Landau gauge, A®'^' = (0, 0, Bxi,0). Note that the counterterms include the 
chemical potential fi and the external field A^^^^. 

To leading order in the coupling constant a = e^/(47r), the fermion self-energy in QED is given by 

S(a:, y) = -Aina-z^Six, y)YD^,{x - y), (2) 

where S{x,y) is the free fermion propagator in magnetic field and D^^{x — y) is the free photon propagator. 

As is well known, the fermion propagator S{x,y) in the presence of an external magnetic field is not translation 
invariant. It can be written, however, in a form of an overall Schwinger phase (breaking the translation invariance) 
and a translation invariant function |25| . i.e., 

S{x,y) =cxp[i<^{x,y)]S{x-y), (3) 

where the Schwinger phase equals $(a;, y) = —eB{xi + yi){x2 — J/2)/2 in the Landau gauge. The Fourier transform of 
S{x — y) is presented in Eq. (jAll) in Appendix \K\ The expression in Eq. ^ implies that the self-energy S(a;,y) has 
an analogous representation 

Y.{x,y)^cyiY>\i<^{x,y)]t{x~y), (4) 

with the same Schwinger phase as in the propagator. 

In this study we use the photon propagator in the Feynman gauge. In momentum space, it reads 

F>p.^{q) = -i'-^ = -i I — 9 ^^ 2 I ■ ^ 9 ^^ K2 I ■ ) ■ (5) 

qi V9o-q -"4 + *^ 95-q2-A2 + iey 

Here we introduced a nonzero photon mass m^ which serves as an infrared regulator at the intermediate stages of 
calculations. Of course, none of the physical observables should depend on this parameter (see Sec. IIVI below). (Note 



that since the classical paper of Stueckclbcrg [2g, it is well known that, unlike non-Abelian theories, introducing a 
photon mass causes no problems in an Abehan gauge theory, such as QED.) We will see in Sec. |IlT] that the leading 
radiative corrections are logarithmically divergent in the ultraviolet region. As in Ref. [131 , we find that the Feynman 
regularization of the photon propagator ([5]) with ultraviolet regularization parameter A presents the most convenient 
way of regularizing the theory. 

The Fourier transform of the translation invariant function S(.t — y) is given by the following expression: 



(2^ 



tip) = -Aina I j—^r Sikh-'D^^ik - p), (6) 



where S{k) is the Fourier transform of the translation invariant part of the fermion propagator and D^u{q) is the 
photon propagator ([S]). 

To linear order in B, the translation invariant part of the free fermion propagator in the momentum representation 
has the following structure: 

5(fc) = 5(°)(fc) + 5(i)(fc) + --- , (7) 

where S^^' is the free fermion propagator in the absence of magnetic field and S^^' is the linear in magnetic field 
part. Both of them are derived in Appendix |X] by making use of a generalized Schwinger parametrization when the 
chemical potential is nonzero. The final expressions for S"'*^' and S^^' can be also rendered in the following equivalent 
form: 

S''\k) ^ , (fco+M)7°-k-7 + m ^3^ 

(Kq + M + *esign(fco))^ — k^ — rn^ 

and 

-S^'\k) = -^WeB (feo + ^h°-fc37^ + m _. (9) 

[(^0 + /i + »e sign(fco))^ — k^ — -m?] 

The self-energy at zero magnetic field 

S(o)(p) = -41^1 ^7^5^°' (fc)7''^M.b-fc) (10) 

determines the counterterms 82 and dm in Eq. ([1]). To calculate the self-energy PH)) . we will use the generalized 
Schwinger parametrization of the fermion propagator S^^'{k), see Eq. (JA11[) in Appendix \K\ Such a representation 
allows a natural separation of the propagator (as well as the resulting self-energy) into the "vacuum" and "matter" 
parts. The former is very similar to the usual vacuum self-energy in QED in the one-loop approximation. The only 
difference will be the appearance of po + M instead of po- The matter part is an additional contribution that comes 
from the (5-function contribution in Eq. (|A11[) . Unlike the vacuum part, the matter one has no ultraviolet divergences 
and vanishes when |/i| < m. 

The explicit expression for the vacuum part reads 



Si°l(p) = £ ^ dx{2m~x [bo + lAf - p • 7] } In ^ 



xK^ 



— x)-m? + x-mi — x{l — x) \(jpQ + /i)^ — p^] 



(11) 



Note that, while the integral over x can be easily calculated, we keep the result in this more compact form. We see 
that the self-energy (jlip becomes identical with the well-known vacuum self-energy in QED in the Feynman gauge 
after performing the substitution po + /i — > po |24l |. Further, using Eq. ()lip . we find that the counterterms in ([T|) are 
defined as follows p4| : 

dp f=m 2tt \2 m? m? 4: j 

5m ^ m - iTLa ^ t^°^^{p) = — 771 ( In ^ + - ) , (13) 



where P — (po + M, p)- Note that the fermion wave function renormalization constant is defined as follows: Z2 = 1 + 5; 



2- 



For completeness, let us calculate the additional matter part of the self-energy due to the filled fermion states given 
by 



Tt'J-^ "J (fco-po)2-(k-p)2 

After performing the integration over the energy and spatial angular coordinates, we find 



e(o)^(p) = -^ f^^^^^kdkh / _ 2m \ ^^ (po + M - Vm2 + fc2)2 _ (fc _ |p|)2 



p|[2V VWT^J (po + M-Vm2 + fc2)2-(fc + |p|)2 

fc(p-7) / fc2 + p2 - (po + M - Vm2 + fc2)2 (po+;x- Vm2 + fc2)2_(fc_|p|)2^ 

1 + JT^^ In- ■ , o , ,oxo 7^^-^ }■ (15) 




|p|Vto2 + /c2 \ 4fc|p| (po + M - \/m2 + fc2)2 - (fc + |p 



While the remaining integral over the absolute value of the momentum k can be also performed, the result will take 
a rather complicated form that will not add any clarity. 

The linear in the magnetic field correction to the translation invariant part of the fermion self-energy in a magnetic 
field reads 

f:M(p)^-A^naJ-^rS^'\k)rD^.{p~k). (16) 

This correction, which in particular contains a chiral shift parameter term, will be analyzed in detail elsewhere [27|. 
However, we use this expression for S'-^-'(p) in the derivation of the leading corrections in the axial current in Sec. IIIII 
below. 



III. THE LEADING RADIATIVE CORRECTIONS TO THE AXIAL CURRENT 

The renormalization group invariant axial current density, which is a quantity of the principal interest in the present 
paper, is given by 

(jf)--Z2tr[7VG(x,x)], (17) 

where G{x, y) is the full fermion propagator and Z2 = 1 + (52 is the wave function renormalization constant of the 
fermion propagator, cf. Eq. ([Ij. 

To the first order in the coupling constant a ~ e^/(47r), the propagator reads 

G{x,y) = S{x,y) + i / d*ud vS{x,u)T,{u,v)S{v,y) + i dud vS{x,u)'Sct{u,v)S{v,y), (18) 

where S{x,y) is the free fermion propagator in the magnetic field, T,{u,v) is the one-loop fermion self-energy, and 
Yjct(u,v) is the counterterms contribution to the self-energy. The structure of the counterterms contribution is 
determined by the last two terms in the Lagrangian density ([T]) . 

In this paper, we make use of the weak magnetic field expansion in the calculation of the axial current density. 
Such an expansion is straightforward to obtain from the general expression in Eq. (|17p and the representation (|18p 
for the fermion propagator. For the fermion propagator to linear in B order, we have 

Six,y) = S^'^\x-y)+te f d^z S^'^^x - z)YS^'^Hz - y)A''J'\z). (19) 



Further, by making use of Eq. (|19p . the weak field expansion of the self-energy follows from the definition in Eq. ([2|). 
(Note that the photon propagator is independent of the magnetic field to this order.) Combining all pieces together, 
we can find the complete expression for the leading radiative corrections to the axial current (J17p in the approximation 
linear in the magnetic field. In this framework, the diagrammatical representation for the leading radiative corrections 
to the axial current are shown in Fig. [1] (for simplicity, we do not display the contributions due to counterterms) . 

Instead of using the expansion for the free propagator in Eq. (J19p , we find it much more convenient to utilize the 
Schwinger form of the fermion propagator ([3]), which consists of a simple phase, that breaks the translation invariance, 
and a translation invariant function. Taking into account that the Schwinger phase ^{x, y) is linear in magnetic field. 



FIG. 1. The leading radiative corrections to the axial current in the linear in magnetic field approximation. Solid and wavy 
lines correspond to the fermion and photon propagators, respectively. Double solid lines describe the axial current insertions 
and the external wavy lines attached to the fermion loops indicate the insertions of the external gauge field. 



we arrive at the following alternative form of the weak field expansion of the fermion propagator in the linear in B 
approximation: 



S{x,y) = 5W(x -v) + z<i>(x,2/)5W(x -y) + S^^^x - y), 



(20) 



where S^^'{x — y) and S^^'{x — y) are the zeroth and first order terms in powers of B in the translation invariant 
part of the propagator. [For the explicit forms of their Fourier transforms see Eqs. ([8]) and ([9]) above.] Of course, the 
representations in Eqs. (fT9|) and ([20|) are equivalent. One can check this exphcitly, for example, by making use of the 
Landau gauge for the external field ^™*. 

Furthermore, Eq. ([2|) implies that a similar expansion takes place also for the fermion self-energy 



S(u,iO = E(°)(u-f)+i$(u,'y)S](°)(u-f) + S(i'(u-w). 



(21) 



The Fourier transforms of the self-energies &^^x — y) and E(^)(a: — y) are given by Eqs. (jTO)) and (fT6|) . respectively. 
Omitting the noninteresting zeroth order in B contribution in Eq. (|18p . we arrive at the following linear in B 
contribution to the propagator: 

G^^\x,x) = 5(1' {x,x)+i f d\d\ \s^^\x ~ u)E(°) (u ~ w)5(") {v - x) + 5(°) (x - u)E(°) {u - v)S^^^ (v - x) 
+ i d ud V 



S^°\x-u)t^^\u-v)S'^°\v-x) 



d'^ud\ [$(a;, u) + $(u, v) + $(w, x)] S'^^^ {x - u)f](°) (u ~ v)S^°^ (v - x). 



(22) 



Noting that $(a;,u) -|- ^{u,v) + ^{v,x) ~ —^ [{xi — Ui){v2 — X2) — (wi — xi)(a;2 — M2)] is a translation invariant 
function, it is convenient to switch to the momentum space on the right hand side of Eq. (P^ . The result reads 



G(i)(cc,.t) = 



d^p 



{2^Y 

eB f d^p 

~2~./ 72^ 



~S^'\p) 



d'^p 

(2^ 



5(i)(p)S(o)(p)5(°)(p) + 5(o)(p)S(o)(p)^(i)(p) + 5(°)(p)s(i)(p)5(°np) 
S(0)(p) 



dpi dp2 dp2 dpi 



(23) 



By substituting this into the definition in Eq. (|17p . we obtain the following expression for the axial current density: 

(24) 
where 



Ol)o 



d^p 



(27r)^ 



rtr 



7^^5(1) (p) 



(25) 



is the contribution to the axial current in the free theory and 



01). = f 



(2^)4 



tr 



7^^ 



dS^'' ip) ^ro) (^) dSi^' iP) ,3 ,5 dS^'' (P) ^ro) (^) dSio^ JP) 



dp 



dp- 



dp2 



dpi 



d^p 

{2^Y 



tr 



7^^5(1) (p)S(o) (p)5(°) {p) + 73755(0) (p)S(") (p)5(i) [p) + 7'7'^^°^ (p)S(i) (p)5(") (p) 



ODct 



(26) 

defines the leading radiative corrections to the axial current. The counterterms contribution (j'Dct in Eq. (pS)) contains 
all the contributions with 82 and 5m- Its explicit form will be given in Subsec. IlII CI below. 



It is instructive to start from investigating the structure of Eq. ((26)) in the free theory (i.e., to the zeroth order in 
a). By making use of the exphcit form of S''^^^(fc) in Eq. ^, we straightforwardly derive the foUowing contribution 
to the axial current density: 



(J5)0 



eB sign(/i) 
4^ 



f d^kS{^i^-k^ -m^) = - 






vV^ 



(27) 



which coincides, of course, with the very well known topological contribution [14J. Note that in contrast to the 
approach using the expansion over the Landau levels, where the contribution to (j'Do comes only from the filled LLL 
states, the origin of the same topological contribution in the formalism of weak magnetic fields is quite different. 
As Eq. ([27]) implies, it comes from the Fermi surface and, therefore, provides a dual description of the topological 
contribution in this formalism. (Interestingly, the origin of the topological contribution in the weak field analysis 
above may have some similarities with the Wigner function formalism J28|.) 

By substituting the propagators ([S]) and (jH]) into Eq. (P5|) , we find the following leading radiative corrections to the 
axial current: 



{jD^ = 32TraeB 



d'^pd'^k 1 

(2^)8 {P-K)l 



{kp + y)[{pQ + ^lf +p\ -pl-rn^] - 2(po + ^)(pifci +^2^2) 

(P2_m2)3(/Y2_^2) 



-2 



(Po + ^J■){Plkl +P2k2 + 2fc3P3 + 4m2) - (fcp + /x)[(po + m)^ + pI + m'^] 

(p2 _ ^2)3(7^2 _ ^2) 

(fco + p)[{po + pf -v\+vl + m^] - 2(]5o + ^^)P3h 



(P2_m2)2(^2_^2) 



{3. 



5/ct 



327raeB 



d^pd^k 1 
(2^)8 {P-K)l 



(fcp + m)[3(po + ^JLf + p2 + m2] - 4(po + A^)(p ■ k + 2m2) 
(P2 _ m?f{K^ - m2) 



(fcp + /x)[3(po + /i)2 - p2 + 3m2] - 2(po + /i)(p ■ k) 

3(p2_ ^2)2(^2 _„,2)2 



ODct. 



(28) 



Here we use the shorthand notation K^ ^ [fcp + /.t + ie sign(fcp)]2 — k2 and P^ = [pp + p + ie sign(pp)]2 — p2. As for the 
definition of {P — K)\, it follows Eq. ([5]). Furthermore, the following replacements have been made in the integrand: 
p\ — > |p2, p2 _j. ip2^ g^j^(;j p^k^ — > i(p • k). These replacements are allowed by the rotational symmetry of the other 
parts of the integrand. 



A. Integration by parts 

It is convenient to represent Eq. (|28| as follows: 



(j|)a = 32^aeB 



d'^pd'^k 



(27r)8 [P-K)] 



4(po + M)[(fco + m)(pp + ju) - P ■ k - 2m2 

(F2_„^2)3(J^2_,^2) 



(fco + m) 



(P2_TO2)2(/i2_^2) 



(fcp + p)[i{po + jif - p2 + 3m2 - 2(p ■ k)] 

3(p2_ ^2)2(^-2 _ ^2)2 



+ 05^>ct. 



(29) 



Since the denominators of the integrand in this expression contain the factors (p2 _ m2)" and (A'2 — rn^)", with 
n = 2, 3, which vanish on the Fermi surface, the integrand in ([29]) is singular there. Therefore, one should carefully 
treat the singularities in the calculation of the axial current. For this, we find it very convenient to use the following 
identity valid for all integer n > 1: 



1 



[[fcp + 11 + ie sign(fcp)]2 — k2 — to2]" [(fcp + /x)2 — k2 — m'^ + ie]" 

-0(|A^|-|fc-o|)0(-fcoM)^("-i)[(fcp+/i)= 



27ri(-l)"-i 



(n-1)! 



c2 - m^l 



(30) 



which can be obtained from Eq. (jA5[) in Appendix 1X1 bv differentiating it n — 1 times with respect to m^. Since the 
first term on the right-hand side has the pole prescription as in the theory without the filled fermion states, we call 



it the "vacuum" part. The second term in this expression takes care of the fiUed fermion states, and we cah it the 
"matter" part. 
One can also obtain another useful relation by differentiating Eq. pop with respect to energy ko, 

d ( 1 \ 2n{kQ + fi) 



dko \ [[ko +^i + ie sign(fco)]2 - m^ - k2]" / [[ko + fi + ie sign(fco)]2 - m^ - k2]"+^ 

^2n^{_^r^Mgin-l) [(fco + ^^r - k^ - m'] [S{ko) S{ko + /x)] ,(31) 

where we made use of Eq. (|30p the second time, albeit with n — ^ n + 1, in order to render the result of derivation 
on the right hand side in the form of the (n + l)th order pole with the conventional ie-prescription at nonzero /j,. In 
addition, we used the following easy to derive result: 

^ [e{\fi\ ~ Ifcol)e(-fcoM)] - sign(M) [<5(fco + /i) - b(ko)\ . (32) 

dko 

We note that 8(ko + /^) in the last term on the right hand side of Eq. ([31]) never contributes. Indeed, this (5-function 
is nonvanishing only when fcg + /i = 0. It multiplies, however, another (5-function, which is nonvanishing only when 
(fco + /i)^ — k^ — 771^ = 0. Since the two conditions cannot be simultaneously satisfied, the corresponding contribution 
is trivial. After taking this into account, we finally obtain 

d ( 1 \ 2n(fco + M) 



dko \ [[fco + M + ^e sign(fco)]2 - vi? - k2]" / [[/j^ + ^^ + jg sign(fco)]2 - m^ - k2]"+^ 

[n — 1)! 

Now, by making use of the above identities, we can proceed to the calculation of {j\)a in Eq. (P^ . We start by 
simplifying the corresponding expression using integrations by parts. Note that the A-regulated representation has 
nice convergence properties in the ultraviolet and, therefore, all integrations by parts in the analysis that follows will 
be perfectly justified. 

The first term in (jDq in Eq. (|29l) is proportional to po + M and contains (P^ — m?)^ in the denominator. This 
suggests to use identity (f33| with n ~ 2 and k ^ p, i.e.. 



jl^ - -i; [jpT^] + '"'■ V - "' - P^ *(P.). m 



Using it, we rewrite the first term in the integrand of Eq. (j29p as follows: 
df^pd^k {ko+ ij){po + ^^|') — p ■'k — 2m? d f — 1 



1st = /i + 32aeS 
= /i + 32aeS 
= /i + 32aeB 



(27r)8 (P - K)\{K^ - m2) dpo \ (P^ - m2)2 

d'^pd'^k 1 d / (fco + m)(po + m) ~ P ■ k — 2m2 



(27r)8 (P2 _ ^2)2 Qp^ y (p _ K)\{K^ - m2) 



d^pd'^k 1 / {ko+^J) (fco+M)(po + M)-p-k-2TO2 d 1 



(27r)8 (P2_to2)2 V(P-ii:)2^(A'2-m2) (A'2-m2) dpo{P-K)l, 

(35) 



where the singular "matter" term, containing the derivative of a (5-function at the Fermi surface, was separated into 
a new function, 

f RA- -i- u f d'^P'^'^k (^0 -h Ai)(po + /^) - P • k - 2m2 . ^ ^ 21, xi \ t'iR\ 

h = 64^- aeBJ ^^-^ i^P - K)l{K^ - m-) ^ [/^ " - " P ] ^(^0). (36) 



We note that the first term in the parenthesis in Eq. (|35l) cancels with the second term in the integrand of Eq. (|29p . 
Then using 

' ' ' ' (37) 



dpo (P - K)l dko (P - K)l 



and integrating by parts, wc find that the sum of first and second terms in the integrand of Eq. (|29p is equal to 

d'^pd'^k 1 f Po + f^ , [n , \/ , \ , o 21 ^ 1 



(1st + 2nd) = h + maeBJ ^^'^^^^ ^^ ^ ^^^^^^^ ^^^^ ( (/°^^2) + [(^o + ^)(Po + m) - P ' k - 2m2] 



d'^pd'^k 1 / (po + m) 0/7 N (^0 +m)(Po + A^) - P ■ k- 2m2\ 



- /i + /2 + 32aeS / ^^^^^ ^^ _ ^^^ (^^^^ _ ^2)2(^2 _ „2) 2(fco + m) ^^2 _ ^2)2(^2 _ ^2)2 



Note that here we used the identity 

d f I \ -2(/co+m) 



dko \ A'2 - m2 / (7^ 



\2 



-2i7r(5(/x2-m2-k2)5(fco), (39) 



which fohows from Eq. ([33]) with n = 1, and introduced another function, which contains the leftover contribution 
with the (5-function, 

f r/1- 2 o f d'^pd'^k{ko + n)ipo + fi)-p-'k-2m'^ , 3 2 1 2^ r/, \ fAn\ 

h = -64*^ aeB j ^^-^ (P _ K)l{P^ - m^f ^^' " "" " ^ ^ °^- ^ ^ 

It is convenient to make the change of variables p ^ k and fc — > p in the first term in Eq. psp . Then, the two terms 
in the integrand can be combined, resulting in 

(1st + 2nd) = .h + h + 32aei?y ^^-^ (P - /v)i(p2 _ ^2)2(^,2 _ „,2)2 • («) 

Finally, by combining the result in Eq. (pij) with the last term in the integrand of Eq. (|29p , we obtain 



64 /■ riV'^fc (^0 + m) 3(P2 - 7^2) + 4p • (p - k) 

3 ' J (27r)8 iP~K)l {P'^-m'^yiK^-m^y 



(Js)" = /i + /2 - -TT-TTaeB / ,- ,g 7^ T7y2—r52 — _2^2^^.^2 — 13:2^2- + 05)ct- (42) 



Using the identity in Eq. ([39]) once again, we rewrite the last expression as follows: 

{jDc = h + /2 + h + {jDct 

64 f dV^fc (fcp-po) / 3 4p.(p-k) 

3 7 (27r)8 (^-^)a V(-P^-™^)(^^-™^) (^^-™^)^(^^-"^^ 



where 

R/l,--,^2„,^R /■ ^4^^4i, q/p2 ^2^ , ^^/^ 1,^ 

.2 _2 1,2 



64i7r2aeB /" d^^pd^k 3(p2 _ ^2) + 4p • (p - k) 

Fi2 



(27r)8 {P - K)l{P^ - m^y 



5{ti^-m'-k^)S{ko). (44) 

Since the first term of the integrand in the second line of Eq. (|43p is odd under the exchange p -^ k, its contribution 
vanishes, and we obtain 

/•3\ f , f , f , /-sn , 64 f d'^pd'^k {ko-pa) 4p • (p - k) 

Finally, by making use of the identity 

and integrating by parts, we derive 

•3x _,,,,,, /-ax ,64 ^ fd^pd^k 2iko-po) ^ -1 



O's)- = /i + /2 + /a + Ot)ct + -TTaeB / ^fr^yrj^ HW^T 5TP ' ^ 



3 7 (27r)8 (P2_ ^2)2(^2 _ ^2)^ ''(P-X) 

f^f^f^/-3\ _^64 ^ /- dV^fc 2(fco-po) 

/i + ,/2 + ./a + 05>ct + —TTaeB / ,„ ,„ — ,,,,,„„ ^r^p • V 



3\ _u^ f dV^fc 2(fco-po) 1 

3™"" i (27r)8 (P-/V)2^(P2_^2)2P' •'(ii:2_^2) 

3, , 64 ^ /"dV^fc 4(fco-Po)p-k 



= ,/i + ,/2 + ./a + (J5>ct + -jnaeB j ^^^^^ ^^ ^^^^^^ _ ^^^^^^^ ^^^^^ 

= h + h + h + ijDcu (47) 



(38) 



where the last integral term in the penultimate line of Eq. (|47|) vanishes because it is odd under the exchange p -^ k. 
Collecting together all contributions, i.e., /i in Eq. ([M]) . /2 in Eq. (|40p and f^ in Eq. (|44p . we have the following 
leading radiative corrections to the axial current: 



{jl)a=6^in^aeB f 



(27r)8 



(fcp + ^l){pQ + /i) - p ■ k - 2to^ p 2 2 21 r/ ^ 
(p,^)^(j^2,„,2) ^ [m -- -P]%0) 



3(po + Ai) - 3(A:o + Af)bo + m) + P - P • k + 3m% . 2 2 



3(P-/Oi(P2-m2)2 



5 (^^ - m^ - k^) 5{h 



05>ct, (48) 



where the first term in the integrand comes from /i , while the second term comes from the sum /2 + /s . The result in 
Eq. (|48p is quite remarkable for several reasons. From a technical viewpoint, it reveals that the integration by parts 
allowed us to reduce the original two-loop expression in Eq. (j29p down to a much simpler one-loop form. Indeed, 
after the integration over one of the momenta in Eq. (j48|) is performed using the (5-functions in the integrand, the 
expression will have an explicit one-loop form. Such a simplification will turn out to be extremely valuable, allowing 
us to obtain an analytic result for the leading radiative corrections to the axial current. 

In addition, the result in Eq. ([iS]) reveals important physics details about the origin of the radiative corrections to 
the axial current. It shows that all nonzero corrections come from the regions of the phase space, where either p or 
k momentum is restricted to the Fermi surface. This resembles the origin of the topological contribution in Eq. (|27p . 
In both cases, the presence of the singular "matter" terms in identities like (|M1) and ([55]) was crucial for obtaining a 
nonzero result. Moreover, by tracing back the derivation of the result in Eq. (|48p . we see that all nonsingular terms 
are gone after the integration by parts. This makes us to conclude that the nonzero radiative corrections to the axial 
current are intimately connected with the precise form of the singularities in the fermion propagator at the Fermi 
surface, that separates the filled fermion states with energies less than /x and empty states with larger energies. 

B. Counterterms contribution 

The calculation of the axial current in Eq. (|48p is still technically quite involved. However, it is relatively straight- 
forward to show [see also the derivation of Eq. (|B7p in Appendix |B] that the right-hand side in ([^5)) without the 
counterterm has a logarithmically divergent contribution when A — > cx3, i.e.. 

To cancel this divergence, we should add the contribution due the counterterms in Lagrangian ([1]). The Fourier 
transforms of the translational invariant part of the counterterm contribution to the self-energy reads 

£(? (p) = 52 [(po + m)7" - P • 7] - '^m , (50) 



where 82 was defined in Eq. (|12p . while 5m = Z2mo — ?n ~ mS2 — Sm and Sm was defined in Eq. (fT3 
We find the following leading order contributions to the axial current density due to counterterms 

/.3\ ./.3\ ,. o /■ d-'p S2{po + fi} .. „ f d^p (po+/x)[<52((po+/^)'-p'+m')-2m(5„J 
(j5>ct = -^2(^5)0 - 4ieB / jT^z^ITlyi Z3T2 " 8«eB ' 

= -SieB 



, (27r)4 (P2 _ ^2)2 y (27r)4 (p2 _ ^2)3 

d^V iPo + m) [S2iP^ - m2) + 2m(m,52 - 6^)] 



(27r)4 (F2 _ ^2)3 

d*p po + fJ. o ■ ^ X jc \ D ^ f d,*p po + fi 



'^^^^^' J (27^4 (p2°_ J^2)2 - 8«"i {^^2 - -5™) eB-, 
eB I-- eBm {17162 - Sm) 



9(m2) J (27r)4 (P2 _ ^2)2 



TT 



2 V/^' ' m^S2 + , =^- (51) 



Here we used the same result of integration as in the topological term, see Eq. (P7)) . 
By making use of the explicit form of the counterterms (|12p and ([T3| , we obtain 

2^-^ \ 2 m2 to2 4 / 47r'^^^2 _ ^2 \^2 m^ 4/ 
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For m ^ |/i|, it reduces to 



^ 27r3 \ 2 m^ m^ 4/ 27r3^ ^2 771^ 4/ 

C. The final result 

The complete expression for the leading radiative corrections to the axial current is given by Eq. (|47p . It consists of 
the counterterm contribution, calculated in the previous subsection, and the additional matter contribution /i+/2+/3- 
The latter is calculated in Appendix IB] For m -^ |/x|, it reads 

aeBu, f k 11 \ aeBm'^ / A 1\ 

^^+^^ + ^^ = ^0-2^+12)+^^0^^2^+6)- (^4) 

Note that this expression has the right ultraviolet logarithmic divergencies (when A -^ 00) that will cancel exactly 
with those in the counterterm ([5^ . Combining the two results, we finally obtain the following leading radiative 
corrections to the axial current in the case m <^ \^\: 

.3^ aeB^i ( ^2^i , ,^ "^^ , 4^ aeBrri^ (^ 2^/2^ 11 



^^■^)" - -^ 1^^^ ^ + ^^ ^ + 3 j - ^^ l^^^ ^ - 12 j • (^^) 

As expected, this result is independent of the ultraviolet regulator A. It does contain, however, the dependence on 
the fictitious photon mass m^. This is the only infrared regulator left in our result. Its origin can be easily traced 
back to the infrared singularity of the wave function renormalization Zi in the Feynman gauge used. As we discuss in 
the next section, this singularity is typical for a class of QED observables, obtained by perturbative methods. As we 
will explain below, in the complete physical expression for the axial current, obtained by going beyond the simplest 
double expansion in the coupling constant and magnetic field, the regulator m^ will likely be replaced by a physical 
scale, e.g., such as \eB\ or a\i^ . 

IV. DISCUSSIONS AND CONCLUSIONS 

Our study of the chiral separation effect in dense QED in the limit of a weak magnetic field suggests a conceptually 
new way to interpret and calculate the axial current density even in noninteracting theory. In contrast to the original 
formulation, which suggests that the topological contribution comes exclusively from the LLL filled states [ij, we 
show that the origin of the same contribution in the formalism of weak magnetic fields (|27| is quite different: it comes 
from the whole Fermi surface. Such a dual description of the topological contribution is of interest on its own. 

Our result for the axial current density, obtained perturbatively in the coupling constant and in linear order in the 
external magnetic field, shows that the chiral separation effect in QED has nonvanishing radiative corrections. To 
leading order, these corrections are shown to be directly connected with the Fermi surface singularities in the fermion 
propagator at nonzero density. This interpretation is strongly supported by another observation: had we ignored the 
corresponding singular terms in the fermion propagator, the calculation of the two-loop radiative corrections would 
give a vanishing result. 

The final result for the leading radiative corrections to the axial current density is presented in Eq. (|55|) . This is 
obtained by a direct calculation of all relevant contributions to linear order in a and to linear order in the external 
magnetic field (strictly speaking, linear in eB because the field always couples with the charge). The result in Eq. ([55)1 
is presented in terms of renormalized (physical) parameters. As expected, it is independent of the ultraviolet regulator 
A, used at intermediate stages of calculations. This is a nontrivial statement since the original two-loop expression for 
the leading radiative corrections contains ultraviolet divergencies. In fact, the divergencies are unavoidable because 
the corresponding diagrams contain the insertions of the one-loop self-energy and vertex diagrams, which are known 
to have logarithmic divergencies. However, at the end of the day, all such divergencies are canceled exactly with the 
contributions due to the counterterms. 

Our analysis shows that the matter contribution, /i + /2 + /a, to the axial current density (calculated in the 
Feynman gauge) has no additional singularities. While functions /i and /2 + /a separately do have additional infrared 
singularities, the physically relevant result for the sum /i + /2 + /a is finite, see Appendix |B] for details. As we see 
from Eq. (|55|) . however, the final result depends on the photon mass tti^, which was introduced as the conventional 
infrared regulator. This feature deserves some additional discussion. 
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It is straightforward to trace the origin of the m^ dependence in Eq. ((55|) to the calculation of the well known result 
for the wave function renormalization constant 82, presented in Eq. (J12p . In fact, this infrared problem is common for 
dynamics in external fields in QED (for a thorough discussion, see Sec. 14 in book [231). The most famous example is 
provided by the calculation of the Lamb shift, when an electron is in a Coulomb field. The point is that even for a light 
nucleus with Za <C 1, one cannot consider the Coulomb field as a weak perturbation in deep infrared. The reason is 
that this field essentially changes the dispersion relation for the electron at low energy and momenta. As a result, its 
four-momenta are not on the electron mass shell, where the infrared divergence is generated in the renormalization 
constant Z2. Because of that, this infrared divergence is fictitious. The correct approach is to consider the Coulomb 
interaction perturbatively only at high energies, while to treat it nonperturbatively at low energies. The crucial point 
is matching those two regions that leads to replacing the fictituous parameter rriy by a physical infrared scale. This 
is the main subtlety that makes the calculation of the Lamb shift quite involved |29[ . 

In the case of the Lamb shift, the infrared scale is related to the atomic binding energy, or equivalently the inverse 
Bohr radius. For smaller energies and momenta, the electron wave functions cannot possibly be approximated with 
plane waves, which is the tacit assumption of the weak field approximation. Almost exactly the same line of arguments 
applies in the present problem of QED in an external magnetic field. In particular, the fermion momenta perpendicular 
to the magnetic field cannot be defined with a precision better than y^eSf , or equivalently the inverse magnetic length. 
This implies that the contribution to the axial current, which comes from the low-energy photon exchange between 
the fermion states near the Fermi surface, should be treated nonperturbatively. Just like in the Lamb shift problem 
[29|, we can anticipate that a proper nonperturbative treatment will result in a term proportional to \Tv{\eB\ / mV) ^ 
with a coefficient such as to cancel the m^f dependence in Eq. (j55|) . 

The additional complication in the problem at hand, which is absent in the study of the Lamb shift, is a nonzero 
density of matter. While doing the expansion in a and keeping only the leading order corrections, we ignored all 
screening effects, which formally appear to be of higher order. It is understood, however, that such effects can be very 
important at nonzero density. In particular, they could replace the unphysical infrared regulator m^ with a physical 

screening mass, i.e., the Debye mass -y/a/i or some combination of powers of a, ^ and ■\/|ei3|. This important issue 
will be addressed elsewhere. 

Another natural question to address is the chiral limit, m — > 0. As one can see from Eq. (|55l). the current (j'Da 
is singular in this limit. This point reflects the well known fact that massless QED possesses new types of infrared 
singularities: Beside the well known divergences connected with soft photons, there are also divergences connected 
with the emission and absorption of collinear fermion-antifermion pairs |30l . l3l| . In addition, because of a Gaussian 
infrared fixed point in massless QED, the renormalized electric charge of massless fermions is completely shielded. 
One can show that this property is also intimately related to the collinear infrared divergences |32|. The complete 
screening of the renormalized electric charge makes this theory very different from massive QED. It remains to be 
examined whether there is a sensible way to describe the interactions with external electromagnetic fields in massless 
QED [13. 

In addition to the quantitative study of the nonperturbative low-energy contributions and the effect of screening, 
there remain several other interesting problems to investigate in the future. Here we will mention only the following 
two. (i) It is of special interest to clarify the connection of the nontrivial radiative corrections to the axial current 
density calculated in this paper with the generation of the chiral shift parameter in dense QED. The preliminary 
analysis [23] shows that there is indeed such a connection but it is more complicated than that in the NJL model 
[la . I21I I22I . (ii) The analysis made in the NJL model shows a lot of similarities between the structure of the chiral 
current in the CSE effect [la, llll |22| with that of the electromagnetic current in the CME one [2j|. It will be 
interesting, therefore, to study the induced electromagnetic current in the CME effect in QED with a chiral chemical 
potential /X5, using the method developed in the present paper. 
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Appendix A: Schwinger parametrization for the fermion propagator at B 7^ and /i 7^ 

The proper-time representation for the ferniion propagator in a constant external magnetic field was obtained long 
time ago by Schwinger [25|. A naive generalization of the corresponding representation to the case of a nonzero 
chemical potential (or density) does not work however. This is due to the complications in the definition of the 
causal Feynman propagator in the complex energy plane when fi ^ 0. The correct analytical properties of such 
a propagator, describing particles above Fermi surface propagating forward in time and holes below Fermi surface 
propagating backward in time, are implemented by introducing an appropriate ie-prescription. In particular, one 
replaces fco + M with kg + fj. + iesign(A:o), where e is a vanishingly small positive parameter. For example, in the 
Landau level representation, the Fourier transform of the translation invariant part of the fermion propagator is 
defined as follows: 

^ ' ^ [^" + /-* + **^ sign(fco)]^ - to2 - A:| - 2n|eB| ' ^ ^ 

where the residue at each individual Landau level is determined by 

D„{k) = 2 [(fco + tJ-h° + m- k^-,^] [P^L„ {2kie) ~ P+L„_i (2fci£2)] ^ 4(fc^ . ^^)lI_^ {2kie) , (A2) 

where -L"(x) are associated Laguerre polynominals. 

Let us start by reminding the usual Schwinger's proper-time representation at zero fermion density, i.e., 



[fco-f iesign(fco)]2-7W2 kl - Ml 



/o 



dse*'*('^«--^"+"), (A3) 



where A^^ = m? + fc| + 2n\eB\. It is important to emphasize that the convergence of the integral and, thus, the 
validity of the representation are insured by having the positive parameter e in the exponent. Unfortunately, such a 
representation fails at finite fermion density. Indeed, by taking into account that 

(A4) 



[/co + M + «esign(fco)]^ - Ml (fco -|- /i)^ - Ml + zesign(fco)sign(fco -I- /i) ' 

we see that the sign of the ze-term in the denominator is not fixed any more. The corresponding sign is determined 
by the product of sign(fco) and sign(fco -I- ^) and can change, depending on the values of fco and ^. For example, while 
it is positive for |fco| > |/j.|, it turns negative when |fco| < |/x| and fco/i < 0. This seemingly innocuous property causes 
a serious problem for the integral representation utilized in Eq. (jA3p . The sign changing «e-term in the exponent 
invalidates the representation at least for a range of quasiparticle energies. 

In order to derive a modified proper-time representation for the fermion propagator, we will make use of the following 
identity: 

1 _ 0(|fco| - \^i\) , ^^1^1 _ i^^i^ ( e{ko^^) , o{~k,fi) 



[ka+fi + iesign{ko)]^- Ml {ko + ^i)'^ - Ml + ie ^"^' ' '^ \{ko + fi)^ - Ml + ie {ko + ^J^ - Ml 

^ -e{\fi\-\k„\)e{~kofi)' ^ ^ 



le 



{ko + ti)^ - Ml + le ''^' ' "" ' '"^^ \{ko+tiy-Ml + te {ko + ^i)^ - Ml - ie ^ 

= {ko + f^r^-M^+ie + ^'''^^'^' " \k,\)9(-kofi)S [(fco + /i)2 - Ml] . (A5) 

The first term on the right hand side of Eq. ()A5P has a vacuum-like ie-prescription and, thus, allows a usual proper- 
time representation. The second term is singular and represents the additional "matter" piece, which would be lost in 
the naive proper-time representation. After making use of this identity, we derive the following modified proper-time 
representation for the propagator: 



o^ /.OO 

S{k) = e^'^i^' ^(-l)"i^„(fc) / ds 



?s[(fco+M) — TTi" — A,'3 — 2n|e_B|+ie] 



Tl = 

CJO 
2 «2 



l/i| - \ko\)ei-ko^i)e-''^' ^(-i)"z?„(fc) 



Tl = 



ds p'^-''li^o+tJ,)^ -m^ —kl-2n\eB\+ie] 



^g ^-is[{ko+fJ,)^ ~m^ -kl-2n\eB\-ie] 



(A6) 
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In order to perform the sum over the Landau levels, we use the following result for the infinite sum of the Laguerre 
polynominals: 

Then we obtain 

/•OO 

Jo 

X [l + 7Vtan(seS)] - 0(1^^1 - |fco|)0(-fcoM) 



X <^ / dse'^[(^°+^) -" -'^3+*^l-''=i^" ^Ms\eB\) ^(^^ ^ ^^^0 + „^ _ k • 7 - (fci72 - ^7') tan(seB)] 

X [1 + 7^7^ tan(seB)] 

+ / rfse-'-'[('^«+^)'-™'-'='-"l+''^i^' ta„(s|es|) ^(^^ _^ ^^-j^o + „j _ k • 7 + (fci72 - k2j^) taii{seB)] 



I 



X [l-7i72tan(seB)] y (A8) 

This is a very convenient alternative representation for the fermion propagator in a constant external magnetic when 
/i ^ 0. It allows, in particular, a straightforward derivation of the expansion in powers of the magnetic field. To 
zeroth order in magnetic field, we obtain 

5(°)(fc) = 5W(fc) + 5i°i(fc), (A9) 

where 

5(°)(fc) = / rfse'"[('=o+A')^-™^-k^+»£l [(fcg + ^)^o + ,„ _ k • 7] (AlO) 

./o 
and 

S°l,{k)=-27T9{\^,\~\ko\)e{-ko^i) [(fco + //)7° + m - k • 7] <5 [(fco + m)' - "i" - k'] (AH) 

are the vacuum and matter parts, respectively. After integration of the proper time and making use of the identity 
in Eq. (|A5[) . we find that this is identical to the usual free fermion propagator ^ in the absence of the field. 

Expanding the expression in Eq. (jASp to linear order in magnetic field, we also easily obtain the following linear in 
B correction to the fermion propagator: 



^(i)(fc) ^ 71726^^ / sdse'*[(''^«+'^^'-"'-'^'+"l + 2nrd{\fi\ - \ko\)9{-kofi)S' [{ka + ^if - m^ 



X [(fco + Ai)7° + m- ^^73] . (A12) 



After integration over the proper time and making use of an identity, obtained from Eq. (jASP by differentiating with 
respect to Ai^^ '^^ obtain Eq. ^. 

Appendix B: Calculation of the /i, /2, and /g terms 

In this Appendix, we give the details of the calculation of the radiative corrections to axial current due to the /i, /2, 
and /a terms. We start from the general form of the result in Eq. ([17|) and calculate separately the two contributions, 
/i and /2 + fs- At the end we combine all contributions and calculate the final result for the /i + /2 + /a contribution. 

1. Calculation of /i 

Starting from the definition in Eq. (|36l) . we find it convenient to rewrite the expression for /i in the following 
equivalent form: 

/i = /i(m^)-/i(A), (Bl) 
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where we took into account that the photon propagator is defined by Eq. ([5]), with A playing the role of the ultraviolet 
regulator. As follows from the definition, 



fiim^) 



~64i7r aeB- 



d 



9(mc)2 J (27r)8 [{P - K'p - mMP^ - m"^) 



S [fi^ ~ml~ k^] S{ko) 



IGiiraeB d 
kp dkp 



kp 



p^dpdp^d^ 



Kpo + m) - pkp^ - 2m2 



(27r)5 (p2 _ p2 _ fc2_ _^ 2pkF^ - "4)[(Po + m)^ - p2 - m2] 



(B2) 



where we integrated over the energy fcg? the absolute value of the spacial momentum k, and all angular coordinates 
except for the angle 9kp between k and p. We also introduced the following short-hand notations: kp = ^ [i^ ~ rn^ 
and ^ — cos 6kp- Note that the auxiliary quantities rric and kp should be replaced by the physical fermion mass m 
and the Fermi momentum pp — -^//i^ — jti^, respectively, at the end of the calculation. 

The integral over the energy pg can be calculated, using the following general result for the energy integration: 



[X{pq + fi)fi + Y]dpo _ vr 



9{n - a) [X{b + ^i)n + Y] e{a - ^) [Xa^^ + (a + b)Y] 



[{b + ^)2 - a2 



a [{a + 5)2 - ^2] 



(B3) 



where a = y^p"^ + rn? and b = 



'2 + fc2 - 2pkp£, 



him 



aeB d 



kpp^dpdS, I 0{pp - p) [/i(6 + /_*) - pkp^ - 2m2] 



'" 2TT^kFdkFj b \ (fe + ^j)2-a2 

d{p - Pf) [afi^ - (a + b){pkF^ + 2m^)] \ 



a [(a + &)2 - /^2] J' 

The integral over the angular coordinate ^ can be easily performed, leading to the following result: 



(B4) 



/i(0) 



aeB d 



pdp 



0{pt 



An^kp dkp 

+9{p-pp)\p + kp -\p- kp\- 




kp j^ (Ai + \P 



kp\ + a){fi + p + kp - a) 



2a 



kl , {a+\p-kp\f -n"^ 

— ^ m — tt: TT 



p- 



[jjL + p + kp + a){^ + \p — kp\ - a) j 

(B5) 



^l- 



Here, without loss of generality, we presented the result only for the case of the vanishing photon mass. This is justified 
because, as wc will see below, the limit m-^ — > does not produce any infrared singularities in the final result for /i. If 
needed, an analogous expression for the case of a nonzero photon mass m^ can be readily written down as well. It can 

be obtained from the above result for making the following three replacements: (i) \p — kp\ 



J{p- kpY + 771?^, (ii) 

p + kp -^ \/(P + kpY + "i2 , and (iii) ^2 _ 3^2 _ ^,2^ _^ ^2 _ ^^ _ ^2^ _ ^^2 g^^- ^^^ places in front of the logarithms. 
After calculating the derivative with respect to kp in Eq. (|B5P and then substituting kp — >■ pp, we obtain 



/i(0) 



aeB 

A'k'^Pp 

+0{p~ 



pdp 



0{PF-P)\ 



2m?p 



PF j^ (m+Pf)^ 



{P^ 



Pf) 



^ _ ^ In ^ 



ipp + H){pp - p^) a [h + pfY -{p- 
2'ni?p\ 2m?p 



PF 



p + pp a{a + p){p'^ - p]r) a{p'^-p%)^ 

It is easy to check that the above expression has a logarithmic ultraviolet divergency, i.e., 

2^2 -I- to2 r ^p 



f\>\0)^aeB 



An^y/li^ — m? 



P 



(B6) 



(B7) 



This confirms that an ultraviolet regularization is required in the calculation. As mentioned earlier, we utilize the 
Feynman regularization (|B1[) . which is equivalent to using the photon propagator in Eq. ([5|). This is the same 
regularization, which is commonly used in the calculation of vacuum diagrams in QED, when the regularized expression 
is obtained from the divergent one by subtracting the contribution with a large photon mass A. In the case at hand. 
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therefore, we need the expUcit expression for the function /i(A). The corresponding calculation is tedious, but 
straightforward. The result reads 



/i(A) 



aeB 



pdp 



6{PF - p) 



P + Pf 



P-PF 



v/(F+P^^F+A2 ^{p - ppY + A2 
{PF -p)(2m2 



^2 ^^2 
A2) 



\/(p - PfY 


+ A2 {2p'p 

{PF 


- 2ppp + A2 + 2ii^ 
+ p)(2m2+A2) 


-Pf)' 


^+A2j 




^J{p + PFY 


+ K\2pl 


+ 2pFP + A2 + 2ii^ 


+ PFf 


^ + A2j 




P., (" + 


^{P + PF? 


+ A2 + aj (^/i + ^{p - 


-PF? 


+ A2-( 


i 



+e{p - pf) 



'pf + A2 + a") (ij. + ^J{p + pFf + k^ 

2p 






{2m' 


+ A2)(p- 


Pf) [a 


+ ViP-PF? 


+ A2j 




ayjip- 


'PF? 

(2m2 


+ A^[2p^-2ppF 

+ A^){p + pp) (^a 


+ A2 + 2a^{p-pp? 


+ A2j 




+ ViP + PF? 


+ A2j 





a^/{p + pp? +A'^ ^2^2 + 2ppp + A2 + 2ay^{p + pp? +aA 
2p2 + 2ppF + A2 + 2a^y{p + pp? +A'^^ 



PF 



In- 



2p2 _ 2ppF + A2 + 2a^(p-pf)2+A2^ 
Finally, as follows from the definition in Eq. (JB1[) . the regularized expression of /i reads 

aeB f^" , I m? m? pp ^ (fi + pp + p + a) {^ + pp — p — a) 



(B8) 



/i 



47r'^fci? Jo 
aeB 



pdp\ 



pdp (2 

PF \ 



{PF - P) {PF + fi) ip + pp) {pp + fl) 
P + PF P-PF 



^In^^ 



{^i + pp - p + a) {fi + pp + p - a) 



An^kp 

rr? {a + p - Pp) 



^W+PF?+^ ViP-PF) 

{2m^ + A^){p - pp) 



A2 



V¥- 



Pf) 



A2 



a{p - pp) {a + p) 


rr? {a + p + pp) 


a{p + pF){a+p) 


Pf , P + Pf Pf 



i-Vip 



PF 



A2 



2a^ip-pp)^+A^) 



A2 ( 2p2 - 2ppp 
{2m^ + A^){p + pp) (a+y^{p + pp)^ +aA 



a\/{p + PF) 
2p2 + 2ppf 



p- pp 



In 
a 2p2 



- A2 ( 2|32 + 2ppp + A2 + 2a^/ip + pp)^+A^ 
A2 + 2ayip 



PF) 



A2 



2ppp + A2 + 2a^(p-pF)2+A2 



(B9) 



A>;j 



Note that, in the first integral below the Fermi surface {p < pp), we took the limit A — ?► 00 because it does not cause 
any problem. It is essential, however, to keep A finite in the second integral above the Fermi surface {p > pp). 

A careful analysis of the regularized expression for /i in Eq. (jB9[) reveals a potentially serious problem: both integrals 
below and above the Fermi surface have infrared logarithmic divergencies, coming from from the regions near pp. 
These divergencies cannot be avoided even when the photon mass is introduced as a regulator. (The divergencies do 
happen to vanish in the theory with massless fermions, m = 0, but this is of no importance as we discuss below.) 
Fortunately, as we will see below, the corresponding divergencies exactly cancel similar infrared divergencies in the 
expression for f2 + fs- Therefore, we come to the conclusion that the appearance of infrared divergencies in /i, as 
well as in /2 + /a, is purely accidental and has no implications on physical observables. They can be viewed as a 
consequence of an ambiguous split of the finite expression /i + /2 + /a into two separate contributions. 
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In order to carefully sort out the cancelation of the above mentioned (unphysical) infrared divergencies, it is useful to 
plicitly separate the divergent terms from 
The divergent part of the expression reads 



explicitly separate the divergent terms from regular ones in the corresponding expression for /i = /} ' '^ + /j '^ 



/i 



(IR,div) 



aeBm 

aeBrri^ 
47r3 



kp—ei 



pdp 







1 



{kp - p){kF + ^i) Jkp + e2 
1 



pdp 






kp + ^ \ £2 



In^-^ 



1 



1 



{p-kp){a+p) 2p2 



2fc, 



In- 



2fc, 



/i 



^^In-^] 
TO^ kp + il 



■BIO) 



where we introduced infrared regulators ei and £2 (with ei, £2 -^ 0) that allow to deal with the problem in a rigorous 
way. Notice that, in the second integral we added a simple regular term, whose only purpose is to insure the ultraviolet 
convergence of the whole expression. The remaining regular part of the expression for /i reads 



n 



(IR,reg) 



aeB 



aeB '•°° 



pdp\ - 



+ 



A-K^k 



pdp\2- 

F JkF \ 

{2rr? + K^){p - k 



{p + kp) {kp + 11) 
p + kp 



kp {^i + kp + p + a) {fi + kp - p — a) 
a {p + kp - p + a) (i-j. + kp + p - a) 

n — kp w? 2m? 



^{p + kpf+K^ ^ip-kp)^+A^ 2p2 



[p - kpY + A2 



^{p-kpY+K'^ 2p2 _ 2pkp + A2 + 2a^/{p-'kpY+K' 



(2m2 + A2)(p + kp) (a + ^{p + kpf+K^^ 

a^ip + kpY + A2 (2p2 + 2pkp + A2 + 2a^{p + kpY + k^^ 

kp , p + kp kp ^ 2p2 + 2pkp + A'^ + 2a J (p + kpV + A2' 

H In In , 

a p- kp a 2p2 _ 2pkp + A2 + 2a^/{p - kp)^ + A2 ^ 

Calculating the integrals in the case ?7i <C |/i|, we arrive at the following results 



a{a + p) {p + kp){a + p) 



A»At 



(Bll) 



n 



(IR,div) 



^^ ln-+ln--l 
Stt'^^ \ £1 £2 



•'^ 27r3 ^ 2n A 



aeBm?' A 

r — m — . 

27r3^ 2/^ 



(B12) 
(B13) 



2. Calculation of /2 + /s 



In this subsection, we calculate the expression for the sum /2 + /a starting from the definition in Eqs. (|40|) and ()44p . 
As we we see below, the corresponding expression has no ultraviolet divergencies. Therefore, we could take the limit 
A — ^ C)o in the expression for /2 + /a, i.e., 



/2 + /3 



GAiTT^aeB d 

3 9K0^ 
32iTTaeBkp d 
3 dijua) 



d'^pd^k 3(po + /i)^ - 3^(po + A*) + P^ - P^fC + Sm^ 



(27r)8 [(P - /v )2 - m2] (P2 _ ,„2) 

p'^dpdpcjd^ 4p2 — p/cj;-^ + 3m2 + 3m2 — 3fi{po + n) 



5 {fi"^ -m^- I?) Siko) 



(27r)5 (p2_p2 



kl + 2pkp^-mmpo + ^iy 



p^ 



(B14) 



where we integrated over the energy fco, the absolute value of the spacial momentum k, and all angular coordinates 
except for the angle 6kp between k and p. We also introduced an auxiliary quantity ma, which should be replaced by 
the physical fermion mass m at the end of the calculation. It should be also noted that the some terms independent 
of TTia were dropped in the integrand of the last expression. This is justified because they vanish anyway after the 
derivative with respect to m^ is calculated. 
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In order to calculate the integral over the energy po, we use again the resuh in Eq. (jB3j). Then, we arrive at 

aeB d [ p^dpd^ ( 0{pF - p)[-'ibiJL + Ap^ - pkpS, - ip^- + im^] 
h + h^ 



67r3 dpF J b \ {b + /.i)^ — a^ 

d{p - Pf) [36/^2 + (a + 5)(4p2 _ ^j^^^ _ 3^2^ ^ 3^2^] ^ 
a [{a + 5)2 - ^2] 



(B15) 



where & = \/p^ + ^f^ ^ 2pkF^ + ml and a = ^Jp'^ +11^— p^ . Note that, in this subsection, we distinguish the quantity 

Pf = ■\//i2 — m2 from the physical Fermi momentum /cj? = \/f^^ ~ m? . After the partial derivative with respect to 
Pf is performed, the value of pf will be replaced by fc^- Similarly, the auxiliary quantity a, which is a function of 
Pf, will be replaced by its physical counterpart y/p^ + m"^ . 

The integral over the angular coordinate ^ can be easily performed, leading to the following result: 



h + h 



aeB d 
12'K^kF dpF 



pdp 



■^F - p)\ P + kF - \p - kF\ + -ipi In 



{pL + \p-kF\f~a^ 
[p+p + kFY — a2 



7(fcg. - p|) + 8p2 + 14to2 {^.i+\p-kF\+a){p+p + kF-a) 



2a 



{^ + p + kF + a){p + \p - kF\ - a) 



, \( , 1 I I, I 2/cFP , . , {a+\p-kF\+ ^){a+p + kF- ^) 

+e{p-pF)\p + kF - \p-kF\ + Vln- — ■ — ■ — -- — ^ 

\ a [a+p + kp + li){a+\p- kp]- p,) 



7{kl-pl) + : 



\Am^ 



2a 



■In 



{a+\p-kp\f- ^P 
(a +p + kpY — ^2 



(B16) 



where, without loss of generality, we again took the limit of vanishing photon mass. 

It is straightforward to show that the integral in Eq. ()B16P has infrared logarithmic divergencies, similar to those 
in function /i . Therefore, we follow the same kind of analysis as in the case of /i and extract explicitly the following 
infrared divergent terms: 



j.(IR,div) ^(IR,div) 

1 2 + /a 



aeB d f ^ 
Un^kp dpp J P^^P 


J 3m\ {f,+ \p^kF[f~a^\ 

0{PF - ei - P In -. — ■ -^-^2 2" 

y fi {fi+P + kpf -a^ ) 




„, .(Zrn?^ a-ii + p-kp 3m'^(kF + n)\ 

+e{p-pF -62) — In + 2 

\ a a + fi + p + kp p^ j 


aeBw? 


[ ' U 


fiei 2kp \ 1 £2 kp + fi 
kp{kp + fi) p J ' kp + fi"' kp + fi kp 




47r3 


kp + fi \ 
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Then, the leftover regular part reads 



p(IR,rcg 



p(IR.rcg 



aeB d 
12Tr-^kF dpp 



pdp 



6{pf — p)\ p + kp - \p — kp\ + Afiln 



(fi + \p-kp\ 



{p,+p + kpY — a2 



7(fcF - p|) + 8p2 + I4m2 {fi+\p-kF\ + a){fi + p + kF -a) im'^ {p,+ \p-kF\) 



2a 



■In 



In^ 



{p+p + kp + a){fi+\p— kp\- a) fi {fi+p + kp) 



/ 2kFp (q+ \p-kp\+ p){a+p + kF- fi) 

+0{p-pF)[p + kp~ Ip-kp] +4/iln— — ■ — - — — 

\ a (a + p + kp + p){a+ \p — kp] - fi) 



nkl - pD + 8p2 + 14m2 ^^ (a + \p - kp])^ - fi^ 



3to2 a-fi + p-kp 3m2(fcF+/i) ., 



2a [a + p -\- kpY — p? a a + fi + p + kp p 

After calculating the derivative with respect to pp in the regular piece (note that a — ^Jp^ + ^2 _ p2^ ^ j q^^ q jg 
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function of pp) and substituting pp — > kp afterwards, wc obtain 



p(IR,rcg) 



p(IR.rcg) 



aeB fkj, 
27r3 V3^ 


^ m^ ^ kp ^ m'^{kp + ^i)\ 
fj. i-iin + kp) 2kj, J 


aeB f pdp 
47r3 J a^ 


j 2m'^p p^ kp{kp + n)+p{a-p)\ 

9{kF - p) —- H In - — r 

y^i[kF + ^i) a kp{kp + n)-p{a+p) j 


1 

+e{p - kp) 
\ 


2pkp 8akp m?{a~p) p^ p — kp m? a + p~kp — fi 
3a 3{a + p) akp + pfi a p + kp a a + p + kp+fi 



;bi9) 



As is easy to check, the integrand is oc 1/p'^ when p — >■ cx). It is clear, therefore, that the expression is convergent and 
no additional ultraviolet regularization is needed. Integrating over the momentum, we finally obtain 



,(IR,rcg) ,(IR,rcg) 

J2 + is 



aeBkp 
27r3 
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■In- 
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3/i kpjjL ^{jjL + kp) 

^2 ,,2 I ^2 



Zkp 



1 



^{kp + /i) 



mM^^^U 



^r + m- kp + /i 
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For m ^ \fi\, the final expressions for the infrared divergent and regular contributions simplify as follows 
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3. Collecting all contributions 



As seen from Eq. (|47p. the final expression for the axial current is given in terms of the sum /i + /2 + /s- The 
corresponding function is obtained by collecting all the divergent and regular terms, calculated in the previous two 
subsections of this Appendix. In the case m <C |/i|, in particular, the result reads 



/l + /2 + h 



/l 

aeB^i 
~2^ 



„(IR,div) „(IR,div) 

+ /2 + /2 
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, A 11 

In \ 
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(B23) 



Notice that all infrared regulators (ei and 62), which were introduced in the divergent parts of /i and /2 + /a cancelled 
out. The only regulator in the last expression is the ultraviolet one A. In the final expression for the axial current 
(|55p . this dependence on the ultraviolet regulator cancels out exactly with a similar dependence coming from the 
counterterms contribution in Eq. ([5^ . 
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